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ABSTRACT 

Current simulations of the rate at which stellar-mass compact objects merge 
with supermassive black holes (called extreme mass ratio inspirals, or EMRIs) 
focus on two-body capture by emission of gravitational radiation. The gravi- 
tational wave signal of such events will likely involve a significant eccentricity 
in the sensitivity range of the Laser Interferometer Space Antenna (LISA). We 
show that tidal separation of stellar-mass compact object binaries by supermas- 
sive black holes will instead produce events whose eccentricity is nearly zero in the 
LISA band. Compared to two-body capture events, tidal separations have a high 
cross section and result in orbits that have a large pericenter and small apocenter. 
Therefore, the rate of interactions per binary is high and the resulting systems are 
very unlikely to be perturbed by other stars into nearly radial plunges. Depend- 
ing on the fraction of compact objects that are in binaries within a few parsecs 
of the center, the rate of low-eccentricity LISA events could be comparable to or 
larger than the rate of high-eccentricity events. 

Subject headings: gravitational waves — relativity — (stars:) binaries: general 
— galaxies: nuclei 



Introduction 



Extreme mass ratio inspirals (EMRIs) of stellar-mass compact objects into supermassive 
black holes are key targets for the Laser Interferometer Space Antenna (LISA). From the 
fundamental physics standpoint, these events are expected to provide the best available 
mapping of the spacetime around a rotating black hole (Ryan 1995, 1997; Hughes 2003). 



1 Department of Astronomy, University of Maryland at College Park, College Park, MD 20742-2421; 
millcr@astro.umd.edu, hamilton@astro.umd.edu, vanessa@astro.umd.edu 

2 Department of Physics and Astronomy, Northwestern University, Evanston, IL 60208; 
m-freitag@northwestern.edu 



- 2 - 



Astrophysically, they may well reveal the numbers of supermassive black holes in a mass 
range (~ 10 5 - 1O 7 M ) that is difficult to probe otherwise (e.g., Greene & Ho 2004). 

Current studies of EMRI rates and properties (Hils & Bender 1995; Sigurdsson & Rees 
1997; Miralda-Escude & Gould 2000; Freitag 2001, 2003; Ivanov 2002; Hopman k Alexander 
2005) have focused exclusively on capture of compact objects by emission of gravitational 
radiation during a close pass. That is, a compact object (for example, a 10 M black hole) 
passes close to the central supermassive black hole (SMBH) and emits gravitational waves 
that shrink its orbit significantly. The black hole then continues to orbit, and if its motion 
is not perturbed significantly by interactions with other stars then it eventually spirals into 
the SMBH. When the orbit becomes detectable with LISA, it has a significant eccentricity 
of typically e ~ 0.5 - 0.9 (Freitag 2003; Hopman & Alexander 2005; but see Ivanov 2002). 

Here we consider a different process, in which a stellar-mass binary containing a compact 
object comes close enough to the SMBH that the binary is tidally separated, leaving one 
object bound to the SMBH and the other almost always ejected to infinity at high speed. 
Tidal separation was discussed recently by Pfahl (2005) as a way to fuel intermediate-mass 
black holes, and has been considered as a method to produce high- velocity stars (Hills 1988, 
1991; Yu & Tremaine 2003; Brown et al. 2005) and as a possible way to deposit high-mass 
stars close to Sgr A* (Gould & Quillen 2003). It was also listed by Hils & Bender (1995) 
and Freitag & Benz (2002) as a mechanism to be examined in the EMRI context, but to our 
knowledge has not yet been explored quantitatively. 

The key point about this process is that, unlike in the two-body capture scenario, no 
energy needs to be dissipated in order to have a capture. As a result, capture can occur 
at much larger radii than is possible in the two-body case: for example, a binary with 
a semimajor axis of tenths of an AU can be captured at pericenter distances of tens of 
AU relative to the SMBH, compared with the ~ 0.1 AU that is required for two-body 
capture. In addition, the semimajor axis of the resulting bound object will be modest, 
perhaps tens of times the pericenter distance (Hills 1991; Pfahl 2005). EMRIs formed in 
this way are therefore relatively immune to perturbations of their orbits that could cause 
them to plunge directly into the SMBH (which lowers rates significantly for EMRIs formed 
by two-body capture; see Hils & Bender 1995; Hopman & Alexander 2005). Combined 
with the higher cross section, this suggests that the overall rate of EMRIs could have an 
important contribution from tidal separation of binaries, even if only a few percent of compact 
objects are in binaries. In addition, the high pericenter after capture implies that when tidal 
separation EMRIs are detectable with LISA they will have eccentricities close to zero (there 
might also be independent paths to low eccentricity, such as production of black holes in an 
accretion disk and their subsequent advection to the SMBH; see Levin 2003). 
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In § 2 we discuss this process in more quantitative detail. In § 3 we list some of the 
questions that will have to be answered to get more specific predictions of relative rates, and 
to interpret LISA observations when they arrive. 

2. Tidal Separation and EMRIs 

2.1. Capture processes 

To evaluate the tidal separation scenario, let us first recall the process of two-body 
capture. Suppose that a point mass m orbits a SMBH of mass M ^> m with an orbital speed 
Voo at apocenter (assumed to be at a large distance). Its orbit will be modified significantly 
if, during its motion, it releases > \mv 1 00 of energy in gravitational radiation. From Quinlan 
& Shapiro (1989), this condition implies a pericenter distance 

r P < r PiGW ~ 0.13 AU(m/10 M ) 2/7 (M/1O 6 M ) 5/7 (v oo /6O km s" 1 )" 4 / 7 . (1) 

We have scaled by 60 km s^ 1 because this is roughly the velocity dispersion inferred for a 
galaxy with a central black hole mass of 10 6 M (Merritt & Ferrarese 2001; Tremaine et al. 
2002; Barth et al. 2005). The time required to spiral into the SMBH would then be much 
less than a Hubble time, except that other stars perturb the orbit significantly (see § 2.2). 
The gravitational radius is r g = GM/c 2 « 0.01 AU(M/1O 6 M ). Therefore, 

r P ,cw/r g « 13(m/10 M ) 2 / 7 (M/1O 6 M Q )~ 2 ' 7 {vjm km s^ 1 )^ 7 . (2) 

For comparison, the radius of the innermost stable circular orbit around a nonrotating SMBH 
is 6r g . As another comparison, detection of an EMRI with LISA will be very difficult if the 
gravitational wave frequency is less than few ~ 2 — 3 mHz, because at lower frequencies there 
is strong unresolvable foreground noise due to double white dwarf binaries in our Galaxy 
(Bender & Hils 1997; Nelemans et al. 2001; Farmer & Phinney 2003). For a circular orbit, the 
gravitational wave frequency is double the orbital frequency (Peters & Mathews 1963). At 
2 mHz, then, the radius of a circular orbit is r(2 mHz) ks 10r 9 (M/10 6 M )~ 2/3 . Therefore, 
a stellar-mass compact object needs to go very deep into the potential well of an SMBH to 
be captured or to be observed with LISA. As a consequence, although the orbit circularizes 
due to emission of gravitational radiation (Peters 1964), the eccentricity in the LISA band 
is still e ~ 0.5 - 0.9. 

Now consider tidal separation. Suppose that a binary with a total mass m and semimajor 
axis a moves towards a supermassive black hole of mass M. If the orbit has a pericenter 
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distance less than 

r tidc ^{^M/m) l ^a 

« 7 AU(M/lO 6 M ) 1 /3( m /ioM )- 1 / 3 (a/O.l AU) , U 

then the binary will be separated by the tidal field of the SMBH. Note that the numerical 
factor in the cube root is correct for a prograde binary on a circular orbit around the SMBH; 
it changes to four for weakly hyperbolic prograde orbits and roughly half this for retrograde 
orbits (Hamilton & Burns 1991, 1992). We scale a by 0.1 AU because such a binary is 
tight enough to survive three-body encounters but wide enough to avoid rapid merger by 
gravitational radiation (see § 3 for further discussion). 

For an initially hard circular binary with component masses 10 M & and 10 M & in a 
hyperbolic pass by a 10 6 M Q SMBH, our numerical simulations suggest that the typical 
eccentricity is e ~ 0.98 after capture, consistent with the results of Hills (1991) and Pfahl 
(2005), who focused on tidal separation of main sequence binaries. For an initial binary 
separation of a = 0.1 AU, the typical pericenter distance after capture is a few AU, and the 
typical apocenter distance is a few hundred AU; both are proportional to the semimajor axis 
of the original binary. We also simulated tidal separation of initially hard circular binaries 
with component masses 10 M and 1 M around a 10 6 M SMBH, representing for example 
a binary with a black hole and a white dwarf or a black hole and a neutron star. We find 
that only a small fraction of encounters lead to ejection of both objects or survival of the 
binary, the rest resulting in capture of one object and ejection of the other. In ~ 40% of the 
captures, the 10 M object becomes bound to the SMBH, with an apocenter distance that 
is a factor of a few larger than for the 10 M — 10 M simulations (as is expected given the 
smaller energy transfer from the 1M object; see Pfahl 2005 for an analytic discussion). In 
the remaining ~ 60% of the captures the 1 M & compact object is captured, which also leads 
to an extreme mass ratio inspiral, but a weaker one than the 10 M — 1O 6 M coalescence 
would produce. 

Because the pericenter distance from binary capture is ^> r g , the orbit circularizes 
dramatically by emission of gravitational radiation and typically has e < 0.01 in the LISA 
sensitivity band, in sharp contrast to EMRIs produced by capture of singles. 



2.2. Effects of nuclear stellar dynamics 

The motion of a binary must be close to radial to be captured. For example, a binary 
with semimajor axis a ~ 1 AU could be captured if it passed within ~ 100 AU of the 
SMBH, but this is tiny compared to the distance of a few parsecs from the SMBH where 
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most binaries presumably lie. It is therefore important to map out some of the dynamical 
processes that will affect the injection into these orbits. These are discussed in detail by 
many authors (e.g., Frank & Rees 1976; Lightman & Shapiro 1977; Magorrian & Tremaine 
1999; Syer & Ulmer 1999), so here we simply quote the results. 

A supermassive black hole of mass M will dominate the dynamics out to the "radius of 
influence" 

r infl = GM/al « 1 pc(M/10 6 M )(6O km s~V CT o) 2 , (4) 

where o"o is the velocity dispersion of stars far outside this radius. At radii r > r inr i, a 
constant velocity dispersion implies a stellar mass density p ~ r~ 2 , whereas at r < r- m R the 
density can take a different slope p ~ r~ 7 , for example 7 = 3/2 or 7 = 7/4 (e.g., Bahcall & 
Wolf 1976; Young 1980). 

For r < ri n fl the orbital time is t Q rb = 27r(r 3 /GM) 1 / 2 , whereas for r > r in fl, t or b = 
27r(r/ri n fl)(GM/o"o). The relaxation time is the time required for the velocity of a star 
to change by of order itself (in magnitude or direction), by deflections due to two-body 
encounters. The local relaxation time for a compact object of mass mco interacting with 
stars of average mass (m) is (Spitzer 1987) 

_ 0.339 a ? '{r) 
trlx[r} ~ In A G 2 (m)m co n(r) " (5) 

Here a(r) is the local velocity dispersion (equal to the orbital speed when r < r infl ), n(r) is 
the local number density, and In A ~ 10 is the Coulomb logarithm. Inside r infl the relaxation 
time is roughly constant. 

For a bound object on a very eccentric orbit, e ~ 1, the angular momentum is much less 
than the angular momentum of a circular orbit with the same semimajor axis. Therefore, 
the angular momentum only needs to change slightly to make an order unity difference in 
the orbit. This timescale is tj(r,e) ~ (1 — e)t r i x (r) (e.g., Hopman & Alexander 2005). 

For a given position R and speed V, the loss cone is defined as the set of directions 

— * 

of the velocity V leading to such small pericenter distances that the object of interest is 
removed from the system. In the full loss cone regime, for which tj < t or b (where tj is 
evaluated for the angular momentum corresponding to the loss cone), objects that enter the 
loss cone and are removed are immediately replaced, within an orbital time, by objects that 
are deflected in from other orbits. In this regime, an object that starts down the loss cone 
is likely to be deflected out of the cone during the orbit. In the empty loss cone regime, for 
which tj > t or b, replacement of objects through the loss cone has to occur over a relaxation 
time. If the number of objects per radius (assuming spherical symmetry) is dN/dr and the 
angle subtended by the loss cone at radius r is #Lc( r )? then the approximate capture rates 
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in the full and empty loss cone regimes are (see Syer & Ulmer 1999) 



dN iun /dr 
dN empty /dr 



6 2 LC (r)(dN/dr)/t orb 

(dN/dr)/ [ln(l/^ c ( r ))t rlx ( r) ] . 



(6) 



Far from the SMBH the loss cone is full, whereas close it is empty. For the binaries, the 
full/empty transition radius (which defines the "critical radius") is comparable to the radius 
of influence, whereas for the singles the transition occurs at ~ 10% of the radius of influence. 
The merger rate is dominated by the region near the critical radius (see, e.g., Frank & Rees 
1976). The smaller critical radius for singles partially compensates for the much larger cross 
section of the binaries, and for 7 = 3/2 the net merger rate enhancement turns out to be 
roughly a factor of ten in favor of the binaries. 

As pointed out by Hils & Bender (1995) and analyzed by Hopman & Alexander (2005), 
there is an additional major effect. A single compact object captured by gravitational ra- 
diation emission typically has a very large apocenter distance, often on the order of tenths 
of a parsec or more. As a result, even after it has first been captured, it has a chance to be 
perturbed in the next orbit. Sometimes, a perturbation will cause the orbit to be so close 
to radial that the object plunges straight into the SMBH. Although this does not affect the 
merger rate, such objects do not contribute to the LISA event rate, because they plunge 
before their orbital period has become shorter than 10 3 ~ 4 s. Hopman & Alexander (2005) 
estimate that ~80-90% of the potential EMRI events are lost in this fashion. This effect 
amounts to reducing significantly the volume from which LI&4-detectable EMRIs can orig- 
inate, which therefore decreases the observed rate. Note, however, that mass segregation 
of black holes into a dense subcluster may reduce the impact of this effect (E. S. Phinney, 
personal communication). 

In contrast, inspirals produced by separation of binaries are not susceptible to this effect. 
The reason is that, as discussed in § 2.1, the apocenter distance is usually only tens of times 
the pericenter distance, hence tj ^> t or b- As a result, we expect that any perturbations will 
be gradual, hence a decrease in the pericenter distance will produce greater gravitational 
radiation emission and thus circularization rather than a plunge. 

Processes that enhance angular momentum diffusion, such as resonant relaxation (Rauch 
& Tremaine 1996) and interactions in triaxial potentials (Holley-Bockelmann et al. 2002; 
Poon & Merritt 2002; Merritt & Poon 2004) will tend to push the full loss cone regime to 
smaller radii, which will enhance rates moderately for both binaries and singles. However, 
assessment of the net effects will require detailed calculations (compare Rauch & Ingalls 
(1998), who show that the total rate of stellar tidal disruptions is at most doubled by 
resonant relaxation, because the rate bottleneck is elsewhere). 
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3. Discussion and Conclusions 

In this paper we have focused on interactions that leave a ~ 10 M black hole in orbit 
around a 10 6 M SMBH. We note that the same process will also enhance rates for neutron 
stars and white dwarfs around 10 6 M Q black holes, likely by an even larger factor, because 
we see from Eq. (2) that direct capture of lower-mass objects at a given speed requires an 
even closer passage to the SMBH. There may also be a moderate effect for intermediate- mass 
black holes (IMBHs) of ~ 10 3 M Q in dense stellar clusters (see Miller & Colbert 2004 for an 
overview of the evidence for IMBHs and their association with clusters), as discussed by Pfahl 
(2005). This could lead to detection with LISA of BH-IMBH orbits (see Baumgardt et al. 
2004 for discussion of gravitational radiation from direct capture onto IMBHs in clusters) or 
orbits of disrupted stars around IMBHs (Hopman & Portegies Zwart 2005), although only 
the nearest sources are likely to be seen (Will 2004). In addition, the loss cone formalism 
is unlikely to be directly relevant here because for M < 10 3 M the wandering radius of the 
black hole is comparable to or larger than its radius of influence (see Merritt 2001) and hence 
the interactions need to be treated as independent binary-single encounters (see Gultekin, 
Miller, & Hamilton 2004 for a recent application to intermediate-mass black holes). 

It is not trivial to estimate the absolute rate of EMRI captures by the mechanism we 
describe, because of the number of processes involved. For example, the fraction of black 
holes in binaries is especially important (see Muno et al. 2005 for a recent discussion in 
the context of the Galactic center). This fraction depends on (1) the fraction of binaries 
that survive stellar evolution (Belczynski et al. 2004), (2) the fraction of those binaries that 
survive interactions in the dense stellar environment of the galactic nucleus (Ivanova et al. 
2005), and (3) the fraction of initially solitary black holes that acquire companions by three- 
body interactions before they are captured by the SMBH. An important output from the 
combination of these processes is (4) the distribution of semimajor axes of binaries containing 
black holes, because very tight binaries (a <0.05 AU for a 10 M — 10 M Q binary) could 
merge by gravitational radiation before being separated by the SMBH, and black holes in 
very wide binaries (more than a few AU) could end up after separation with semimajor axes 
so large (more than ~ 0.1 pc) that perturbations during a single dynamical time drop them 
into plunge orbits, making them undetectable with LISA (Hopman & Alexander 2005). 
Mass segregation will tend to move black hole binaries to regions of higher density and 
higher velocity dispersion, where three-body interactions are important, hence we must also 
compute (5) the evolution of the semimajor axis and companion mass as a function of time, 
versus the probability of capture by the SMBH as a function of time, to estimate the true 
distribution of semimajor axes after capture, and hence the subsequent evolution of the 
orbit under the influence of relaxation and gravitational radiation. All of these processes 
will require careful computation in future work. 
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Without knowing the absolute rate, we can parameterize the ratio of the rate of captures 
due to binary separation to the captures due to singles as 

-^binary /b^binary (/binary, LISa) 

-^single f s ^"single (/single.LISA) 

where fb is the fraction of compact objects that are in binaries that are neither too tight 
nor too wide (see above); f s is the fraction that are single; 7^bmary is the total rate of tidal 
separations per binary; 7^ sing i c is the total rate of gravitational radiation captures per single; 
(/binary,LiSA) is the overall fraction of binary sources captured in orbits tight enough to spiral 
into the LISA band within a Hubble time; and (/ S mgic,LisA) is the overall fraction of captured 
singles that end up detectable with LISA (rather than being perturbed into plunge orbits). 
Our current best guesses are 7^ b inar y /Single ~ 10 and (/binar y ,LisA)/(/singic,LisA) ~ 1 - 10. 
Therefore, if the steady-state binary fraction (reduced by merging at small semimajor axes, 
and by ionizations at large semimajor axes) is fb > 0.01 — 0.1, EMRIs from binaries could 
dominate the total rates. 
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